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Abstract. This paper is devoted to the study of the Hausdorff dimension of the singular set of 
the minimum time function T under controllability conditions which do not imply the Lipschitz 
continuity of T. We consider first the case of normal linear control systems with constant coefficients 
in M. N . We characterize points around which T is not Lipschitz as those which can be reached 
from the origin by an optimal trajectory (of the reversed dynamics) with vanishing minimized 
Hamiltonian. Linearity permits an explicit representation of such set, that we call <S. Furthermore, 
we show that S is 'H JV_1 -rectifiable with positive 'H JV_1 -measure. Second, we consider a class 
of control-afnne planar nonlinear systems satisfying a second order controllability condition: we 
characterize the set <S in a neighborhood of the origin in a similar way and prove the 'H 1 -rectifiability 
of S and that K X (<S) > 0. In both cases, T is known to have epigraph with positive reach, hence 
to be a locally BV function (see |13l 115] ). Since the Cantor part of DT must be concentrated in 
<S, our analysis yields that T is SBV, i.e., the Cantor part of DT vanishes. Our results imply also 
that T is locally of class C 1 ' 1 outside a 'H JV_1 -rectifiable set. With small changes, our results are 
valid also in the case of multiple control input. 



1. Introduction 

Consider the control system 

{x = F(x) + G(x)u, Id < 1, x e R N , 
x(o) = e 

The minimum time function T(£) to reach the origin from £ under the above dynamics is well known 
to be, in general, both non-smooth and non-Lipschitz. Several papers were devoted to the partial 
regularity of T. In particular, we quote results devoted to establishing (semi) convexity /concavity 
properties of T under various assumptions (see jSJ El Q21 E51 E21 E]), together with [3J, which is 
concerned with planar systems. 

In this paper we concentrate mainly on the lack of Lipschitz continuity of T, which is essentially 
due to the lack of first order controllability. More precisely, even if at some x the right hand side of 
(|1.1|) does not point towards the origin, i.e., the scalar product between any vector in F(x) + G(x)U 
and x is merely nonnegative, it is still possible that a trajectory through x reach the origin, provided 
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the Lie bracket [F(x),G(x)] has nonvanishing scalar product with the missing direction x (higher 
order controllability). The price to pay is a slower approaching to the origin: one needs to switch 
between G and —G. like a sailor which has to beat to windward. The simple example x = u 6 [—1, 1] 
exhibits this behavior: at every point of the xi-axis (we set x\ = X2, £2 = u) the right hand side of 
(jl.ip is vertical and T is not locally Lipschitz in the whole of W 1 (but at the points of the xi-axis 
is Lipschitz). By introducing the minimized Hamiltonian 

h(x,0 := (F(x),C)+min(G(xKC), 

the condition of non-pointing towards the origin becomes 

h(x,() > 0, 

where C is a normal to the sublevel of T corresponding to T(x). Since the minimized Hamiltonian is 
constant and nonpositive along every optimal trajectory, it is natural to expect that non-Lipschitz 
points of T lie exactly where such Hamiltonian vanishes. In fact, in Section 13.21 we prove this 
characterization. 

Let S be the set of non-Lipschitz points of T. In Sections [Hand Owe characterize S using points 
which belong to an optimal pair (i.e., an optimal trajectory together with a corresponding adjoint 
arc) of (jl.ip with vanishing Hamiltonian, and, for the linear case, we give an explicit representation 
of S. As a consequence, we show that at each x G S the sublevel TZt(x) 1S tangent to S, in the sense 
that there exists a normal vector to TtT(x) a ^ ^ which is tangent to the optimal trajectory reaching 
x from the origin. The result is valid for both normal linear systems with constant coefficients in 
any space dimension (see Theorem I4.4p and for smooth nonlinear two dimensional systems such 
that the origin is an equilibrium point, the linearization at is normal and furthermore DG(0) = 
(see Theorem I5.2p . The condition DG(0) = ensures that the nonlinearity is sufficiently mild to 
preserve a linear like behavior in a neighborhood of the origin whose size can be estimated. In both 
cases it is known that the epigraph of T has locally positive reach (see [13\ LT5] ) . Reasons for the 
restriction to two space dimensions in the nonlinear case are discussed in the paper [15], to which 
the present work owes some results. 

Our first result is the 'H Ar_1 -rectifiability of S for the linear single input case, see Theorem 14.51 
and, respectively, the T^-rectifiability for the nonlinear two dimensional case, see Theorem 15. 21 For 
the linear case, the switching function 

g ( (t) = (t,e At b), 

where b is a column of the matrix B, plays an important role. Actually we partition S according to 
the multiplicity of zeros of and embed each part into a locally Lipschitz graph. The nonlinear case 
is handled by showing that S consists of optimal trajectories with vanishing Hamiltonian. Since, 
due to the space dimension restriction, such trajectories are at most two for the single input case, 
the T^-rectifiability is clear. We observe that the investigation of the regularity of the Minimum 
time front, i.e., the boundary of IZt, performed in [31 Chapter 3] cannot provide information on 
non-Lipschitz points of T, since an analysis of sublevels is not enough to describe the epigraph of 
a function. 

These rectifiability results shade some light on the propagation of singularities for minimum time 
functions. In fact, the positive reach property of epi(T) implies that T is locally semiconvex outside 
the closed set S (see Theorem 5.1]). The structure of singularities of semiconvex functions is 
well understood (see {9l Chapter 4]): in particular, a locally semiconvex function is of class C 1 ' 
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outside a T-L -rectifiable set. Therefore our rectifiability results for S imply that T is of class C > 
outside a closed ^ Ar_1 -rectifiable set. We observe that, for general functions whose epigraph satisfies 
a uniform external sphere condition, the Hausdorff dimension of the set of non-Lipschitz points was 
proved to be less or equal to n — 1/2, with an example showing the sharpness of the estimate (see 
|21|. Theorem 1.3 and Proposition 7.3]). The present paper therefore improves that result, for the 
particular case of a minimum time function. We prove also (see section [7]) a converse (propagation) 
result: for 7~L N ~ 1 -a.e. x G S such that T(x) is small enough there exists a neighborhood V such 
that H N ~ 1 (S n V) > 0. In particular, for single input normal linear systems we show that, for all 
t > small enough, the set S D IZt, up to a 'H Ar_2 -negligible subset, is a C 1 -surface of dimension 
N — 2. Apparently, this is the first result in the literature concerning propagation of non-Lispchitz 
singularities. 

Applications of the above results to time optimal feedbacks will be discussed elsewhere. 

The positive reach property of epi(T) implies also that T has locally bounded variation (see 
[111 Proposition 7.1]). Therefore, as a consequence of the above analysis we obtain that T belongs 
to the smaller class of locally SBV functions, namely the Cantor part D C T of its distributional 
derivative, which is a Radon measure by definition of BV, vanishes. In fact, on one hand D C T must 
be concentrated on the set S of non-Lipschitz points of T, on the other the rectifiability properties 
that we proved for S yield exactly the SBV regularity of T. To our best knowledge this property 
of T is observed here for the first time. 

2. Preliminaries 

2.1. Nonsmooth analysis, sets with positive reach, and geometric measure theory. 

The space dimension is denoted by N and we suppose N > 2. The unit sphere in is denoted 
byS*" 1 . 

Let K C l w be closed with boundary bdry.fr. We need some concepts of nonsmooth analysis (see, 
e.g. 10. Chapters 1 and 2]). Given x G K and v G 1^, we say that v is a proximal normal to K 
at x, and denote this fact by v G Nk(x), provided there exists a = a(v,x) > such that 

(v, y — x) < a \\y — x\\ 2 , for all y G K. 

If K is convex, then Nk(x) coincides with the normal cone of Convex Analysis. The set of limiting 
normals to K at x is denoted by N^(x), and consists of those v G M. N for which there exist 
sequences {xj}, {vi} with xi — > x, Vi — > v, and Vi G Njc(xi). The Clarke normal cone N^-(x) equals 
coN^(x), where "co" means the closed convex hull. 

Let Q C M. N be open and let / : Q, — > K U {+00} be lower semicontinuous. The epigraph of / 
is epi(/) = {(x,y) £(lxl:|/> fi x )}- The proximal subdifferential df(x) of / at a point x of 
dom(/) = {x : f(x) < +00} is the set of vectors v G M. N such that 

(v,-l) G N cpi{f) (x,f(x)). 

The horizon subdifferential d°°f(x) of / at a point x G dom(/) is the set of vectors v G M. N such 
that 

(v,0) G iV epi(/) 0r,/(x)). 

This concept is connected with the lack of Lipschitz continuity of / around x (see, e.g., [24, Chapter 
9]) and will be used mainly in Section 3. 
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Sets with positive reach will play an important role in the sequel. The definition was first given by 
Federer in [IB] and later studied by several authors (see the survey paper [IB]). 

Definition 2.1. Let K C M. N be locally closed. We say that K has locally positive reach provided 
there exists a continuous function (p : K — > [0, +oo) such that the inequality 

(2.1) (v,y - x) < ip(x)\\v\\ \\y - x\\ 2 
holds for all x,y G K and v G Nk{x). 

In particular, every convex set has positive reach: it suffices to take <p = in (|2.ip . 
Continuous functions whose epigraph has locally positive reach will be crucial in our analysis. Such 
functions enjoy several regularity properties, mainly studied in [11] . We list two of them which will 
be used in the sequel. 

Theorem 2.2. Let C M> N be open and let f : ft — >• K be continuous and such that epi(f) has 
locally positive reach. Then 

(i) f is a.e. differentiable in $7, 

(ii) f has locally bounded variation in fL 

A brief survey of basic notions on functions with bounded variation will be presented in Section 
6. Our main reference, also for other basic concepts in geometric measure theory, is [TJ. 

2.2. Control theory. Consider the following autonomous control system 

' y(t) = f(y(t),u(t)) a.e., 

(2.2) I u{t) G U a.e., 

y(0) = x, 

where the control set IA C M. M is nonempty and compact and / : M. N x IA — > M. N is continuous 
and Lipschitz with respect to the state variable x, uniformly with respect to u. We denote by U a d 
the set of admissible controls, i.e., all measurable functions u such that u(s) £ IA for a.e. s. Under 
our assumptions, for any u(-) G U a di there is a unique Caratheodory solution of (|2.2p denoted by 
y x ' u (-). The solution y x ' u (-) is called the trajectory starting from x corresponding to the control 

«(•)■ 

We will focus mainly on control systems which are linear or nonlinear with respect to the space 
variable and linear and symmetric with respect to the control. More precisely, we will consider the 
linear control system 

y(t) = Ay(t) + Bu(t) a.e., 

(2.3) I u{t) G U = [-1, 1] M a.e., 

, y(0) = x, 

where 1 < M < N and A G Mtvx7V> B G MjvxM and IA = [-1, 1] M B (m, . . . ,u M ) =■ u, together 
with the nonlinear two dimensional control system 

' y(t) = F(y(t)) + G(y(t))u(t) a.e., 

(2.4) I u{t) G [-1,1] M a.e., 

k y(0) = x, 

where F and G are suitable vector fields (the actual assumptions will be stated later) and 1 < 
M < 2. We will use also the notation B = (pi, . . . , 6a/) or G = (gi, . . . , <7m)j where each entry is 



NON-LIPSCHITZ POINTS AND THE SBV REGULARITY OF THE MINIMUM TIME FUNCTION 



5 



an iV-dimensional column. Note that x is reachable by a solution of (|2.3p at time t if and only if 
the following (equivalent) conditions hold: 

(2.5) x = e At x + [ e A{t - s ^ Bu(s) ds and x = e~ At x - I e- As Bu(s)ds, 

Jo Jo 

where u(-) £ U a d- 

For a fixed x £ M. N , we define 

9{x,u) := min {t > | y x ' u (t) = 0}. 

Of course, 9{x,u) £ [0, +oo], and 9(x,u) is the time taken for the trajectory y x ' u {-) to reach 

0, provided 9(x,u) < +oo. The minimum time T(x) to reach from x is defined by T(x) := 
inf {9(x, u) | u(-) £ ^/ a d} an d under standard assumptions the infimum is attained. A minimizing 
control, say u(-), is called an optimal control. The trajectory y x ' u (-) corresponding to u(-) is called 
an optimal trajectory. 

Denote by TZt the set of points which can be steered to the origin with the control dynamics 
(|2.2p within the time t. Then TZt is the set of points which can be reached from the origin with the 

reversed dynamics 

x(t) = -f(x(t),u(t)) a.e., 

(2.6) < u(t) G U a.e., 

k i(0) = 0, 

1. e., TZt is the sublevel {x G 1^ : T{x) < i) of T(-). If u is an admissible control steering x to the 
origin in the minimum time T(x), then the Dynamic Programming Principle implies that for all 
< t < T(x) the point y x,u (t) belongs to the boundary of TZt- 

We state now Pontryagin's Maximum Principle, giving first its linear version for the special case 
we are interested in. 

Theorem 2.3. Consider the problem (|2.3p under the following (normality) condition: 
for every column hi of B, i = 1, M , we have 

iaDk[b i ,Ab i ,...,A rf - 1 bi] = N. 

Let T > and let x 6 M. N . The following statements are equivalent: 

(i) x £ bdryT&T; 

(ii) there exists an optimal control u steering x to the origin in time T; in particular, T{x) = T; 

(iii) (Pontryagin Maximum Principle) for every ( G Nn T (x),C ^ 0, we have 

(2.7) Ui (t) = -sign ((C,e~ A %)) a.e. t £ [0,T], 
for all i = 1,2,..., M. 

A well known reference for this result is |20i Sections 13 - 15]. 

Remark 2.4. If the system A2.3\) is normal then (A, B) satisfies the Kalman rank condition. There- 
fore the minimum time function is everywhere finite and continuous (actually Holder continuous 
with exponent 1/N, see, e.g., Theorem 17.3 in |2L)j and Theorem 1.9, Chapter IV, in [2] and refer- 
ences therein). 
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Before stating Pontryagin's Principle for the nonlinear case (|2.2p . we need to introduce the mini- 
mized Hamiltonian. We define for every triple (x,p,u) el^xK^x [-1, 1] M , the Hamiltonian: 

H(x,p,u) = {p,f(x,u)) 

and the minimized Hamiltonian: 

h(x,p) = mm{H(x,p, u) : u G [— 1, 1] M }. 

Observe that if x is steered to the origin with respect to the system (|2.2p by the control u(-) in 
the time T, then the origin is steered to x with respect to the reversed dynamics (|2.6[) in the 
same time T by the control u(t) := u(T — t). The corresponding trajectory will be denoted by 
y(t) := y x ' u (T — t). Then Pontryagin's Principle reads as follows. 

Theorem 2.5 (Pontryagin's Maximum Principle for nonlinear systems). Fix T > and let x G 

together with an optimal control steering x to the origin in the time T. Then there exists an 
absolutely continuous function A : [0,T] — > M. N , never vanishing, such that 

(i) \{t) = \(t)D x f(y(t),u(T-t)) } a.e. t € [0,T\, 

(ii) H (y(t), X(t),u(T -t)) = h (y(t), X(t)), a.e. t G [0, T], 

(iii) h (y(t), A(i)) = constant, for all t E [0, T], 

(iv) \(T)£NC t (x). 

This formulation of the Maximum principle can be obtained by using the classical one (see, e.g., 
[251 Theorem 8.7.1]) for the reversed dynamics. Observe that the condition (|2.7|) is equivalent to 
the minimization condition (ii) for the case of linear systems. 

3. Properties connected with the minimized Hamiltonian 
This section is mainly technical and consists of two subsections. 

3.1. Minimized Hamiltonian and normals to epi(T). This section is concerned with a relation 
between normals to the sublevels of T and normals to epi(T), which was one of the main tools used 
in [131 02] in order to prove that epi(T) has positive reach. We give here a unified and slightly 
generalized presentation, in order to use it in the sequel. 

We recall first that a point x G W N \ {0} is defined to be an optimal point for (|2.4p if there exist 
x\ such that T(xi) > T(x) and a control u with the property that y Xl ' u [-) steers x\ to the origin in 
the optimal time T(x\) and y Xl > u (T(xi) — T(x)) = x, i.e., if there exists an optimal trajectory for 
(|2.4p which passes through x. It is easy to see, using Pontryagin's Maximum Principle, that for a 
normal linear system every point is optimal. Indeed, it is enough to extend the adjoint vector and 
choose a control which maximizes the Hamiltonian (for the reversed dynamics). For nonlinear two 
dimensional systems, sufficient conditions for this property were given in |15^ Theorem 6.2]. 
Second, we recall that sublevels of T are convex in the case of linear systems (see, e.g, |20[ Lemma 
12.1]) and so, in particular, they have positive reach. For nonlinear two dimensional systems, 
sufficient conditions for this property were given in [151 Theorem 5.1(c)]. 

After the above preliminaries, we can state our result. The assumptions are indeed strong, but 
we emphasize the fact that they are all satisfied in the two cases we are going to consider in the 
present paper (see Theorem 3.7 in [13] and Theorem 6.5 in |15j). 

Proposition 3.1. Consider the general control system \2.2t) with the following assumptions: 
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(i) hi C M M is compact and {/(x, u) : u G is convex for every x G M. N . 

(ii) / : l w X W — )■ is continuous and satisfies 

\\f{x,u)-f(y,u)\\<L\\x-y\\, Vx, y G M", u G W, 

/or a positive constant L. Moreover, the differential of f with respect to the x variable, 
D x f, exists everywhere, is continuous with respect to both x and u and satisfies 

\\D x f(x,u) - D x f(y,u)\\ < Li\\x - Vx,y G R^^ G U, 

for a positive constant L\ . 

Let x G R N \ {0} and let T(x) be the minimum time to reach the origin from x. Assume that there 
exists a neighborhood V of x such that 

(1) T is finite and continuous in V, 

(2) every y G V is an optimal point, 

(3) for every y G V the optimal control steering y to the origin is unique and bang-bang with 
finitely many switchings, 

(4) there exists r > T(x) such that !Z t has positive reach for all t < r. 

Let C G Nn T{x) {x). Then 

(a) h(x,C) < 0. 

(b) (C,h(x,0) GiV cpi(r) (x,T(x)). 

Proof. Since £ G Nn T{x) (x) and TZ>t(x) nas positive reach, there exists a constant a > such that 

(C,y-x) <a\\C\\\\y-x\\ 2 , 

for all y G K T (x)- 

Let (x(-),n(-)) be an optimal pair for x and let A : [0,T(x)] — > R N be absolutely continuous and 
such that 



(3.1) 



I A(i) = A(t)^/(x(T(x)-t),w(r(x)-t)), 
1 A(T(x)) = C- 

and h(x(t),X(t)) = (X(t), f(x(T(x) - t),u{T(x) - t))) = h(x,() for all t G [0,T(x)] (see Theorem 
3.1 in [3]). 

Observing that x(t) G 7Zt(x) for all < i < T(x), we have 

(Cx(t)-x) <o-||C||||x(t)-x|| 2 . 

By using GronwalPs lemma, there is a suitable Mi > such that \\x(t) — x|| < M\t. Since / is 
Lipschitz with respect to x, we have for some M2 > and for all < t < T(x), 

(C,j J' f{xMs)ds)<M\\C\\t. 
Taking the upper limit, we obtain 

limsup/e, - / f(x,u(s)ds\<0, 

which implies 

<C,/(x,n)> <0, 

for some u G ZY, since {/(x, u) : u G li\ is convex and U is compact, whence h(x, Q < 0. 
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We are now going to show that (£, h(x, ()) G N ep UT\(x, T(x)) : i.e., there is a a > such that 
(3.2) ((C,0), (y,/3) - (x,T(x)) < a(\\y - x\\ 2 + \T(x) - /3| 2 ), 

for all (y,/3) in a neighborhood of (x,T(x)), say G V x [0, r], where i? := h(x,Q), /3 > T(y), 

and r > T(x) is such that T(y) < r for all y G V. There are two possible cases: 

(i) T(y) < T(x), 

(ii) T(y) > T(x). 

In the first case, since y G 1Zt(x) an d T^t{x) nas positive reach, there is a i^i > such that 

((,y-x) ^^HCIIIIy-xll 2 . 

Since ■d = h(x,C) < 0, if /3 > T(x) then (|3.2p is satisfied. If instead /3 < T(x), then we set 
x% = x(T(x) — (3). By Gronwall's lemma, there is some K > such that 

\\x -xi\\ = \\x - x{T(x) -P)\\< K\T{x) - p\. 

We have 

((,y-x) = (A( i 9) jI /-xi) + (A(T(x))-A(y3),y-xi) + (C,xi-x) 
= :(/) + (//) + (///). 

We now consider (/). Since y G Hp and A(/3) G (x\) (see Definition 2.3 and Corollary 4.8 in 
[19]), owing to the fact that Hp has positive reach, there exist K2,K^ > such that 



(/) < ^2||AC8)||||y-xi 



|2 



< 2K 2 ||A(/3)||(||y-x|| 2 + ||x-x 1 || 2 ) 

< K 3 (\\y-xf + \T(x)-/3\ 2 ). 

Let us now consider (//). We have, for suitable constants K^,K^ > 0, 

(//) < ||A(r(x))-A(/3)|||| 2/ -x 1 || 

< K4T{x)-0\(\\v-x\\+K\T(x)-P\) 

< K B (\\y-x\\ 2 + \T(x)-P\ 2 ). 

Finally, we have, for a suitable constant Kq > 0, 

(III) = F* \x(T(x))J(x(s),u(s))ds 
Jo 

T(x)-/3 rT(x)-P 

(X(T(x) - s), f(x(s),u(s))ds + / <A(T(x)) - X(T(x) - s), f(x(s),u(s))ds 



(since the minimized Hamiltonian is constant) 

f T(x) 

< (T(x)-P)h(x,C)+ / K e \T(x) - s\ds 

h 

= (T(x)-P)h(*,0 + !y\T(x)-f3\ 2 . 
Putting the estimates together, we obtain 

(C, y-x)< (T(x) - P)h{x, C) + K 7 (\\y - x\\ 2 + \T(x) - /3| 2 ), 
for a suitable positive constant Kj. The proof of (|3.2p is concluded in the case (i). 
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We are now going to consider the case (ii). Since i? < 0, it is enough to prove (|3.2p for (3 = T(y). 
Since x is an optimal point, there exists x\ such that T{x\) = T(y) together with an optimal pair, 
still denoted x(-) and u(-), such that x{T(y) — T(x)) = x. Let A(-) denote the extension up to the 
time T(y) of the solution of (|3.ip . Since the optimal control is unique and bang-bang with finitely 
many switchings, it is easy to prove that h{x{T(y) — t), X(t)) = (X(t), f(x(T(y) — t),u(T(y) —t))) = 
constant for all t G [0,T(y)]. Then by using the same argument of the case (i), one can easily show 
that (|3.2p holds true. The proof is complete. □ 

3.2. Minimized Hamiltonian and non-Lipschitz points. This section is devoted to identify 
points around which the minimum time function T is not Lipschitz as points where the proximal 
normal cone to epi(T) contains a horizontal vector £ ^ 0. It will also turn out that if x is a non- 
Lipschitz point and £ is such a vector, then h(x, £) = 0. A kind of converse statement can also be 
proved. 

All results are valid in the domain where epi(T) has positive reach. So, for linear systems, 
they hold globally, while for nonlinear two dimensional systems of the type (|2.4p . they hold in a 
neighborhood of the origin whose size can be estimated and depends only on the data. 

Definition 3.2. We say that a function T : R-^ — > R is non-Lipschitz at x provided there exist two 
sequences {xi},{yi} such that Xi / yi for all i, {xi},{yi} converge to x and 

y \T( yi )-T(xj)\ 

Inn sup 1| r. = +00. 

Observe that the set of non-Lipschitz points is closed. 
The first result does not require T to be a minimum time function. 

Proposition 3.3. Let C R w be open and let T be continuous in Q, and such that epi(T) has 
locally positive reach. Let x G £1. Then T is non-Lipschitz at x if and only if there exists a nonzero 
vector ( £ R^ such that 

(CO) £ N epi{T) (x,T(x)). 

Proof. By Theorem 9.13 in [23], T is non-Lipschitz at x if and only if d°°T(x) contains a nonzero 
vector Q. This condition is equivalent to (CO) G N epi ^(x,T(x)). □ 

Now we restrict ourselves to the case where T is the minimum time function to reach the origin 
for (|2.3p or for (|2.4p . We assume that the conditions ensuring that epi(T) has positive reach are 
satisfied (see [131 Theorem 3.7] and [151 Theorem 6.5]). 

Proposition 3.4. Let T denote the minimum time function to reach the origin for \2. 3\) or for 
pO ). Let and 5 > be such that the epigraph of T restricted to B(x,S) has positive reach. 

Let C G 1^, C + 0. Then 

C G d°°T(x) if and only if h(x, () = and ( G A r ^ T(s) (x). 

Proof. Recalling Proposition 13.31 ( G d°°T(x) if and only if (CO) G N cpi ^(x,T(x)), i.e., for a 
suitable constant c > 0, 

(3-3) (C y-x)< clldl {\\y - *|| 2 + |/3 - T(x)\ 2 ) , 

for all y G B(x, 5) and for all (3 > T{y). 
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Let C G d°°T(x). If y G 71t(x)i then we can take = T(x) in fj3.3[) and so 

(C,y-^)<c||C||-||y-s|| 2 , 

i.e., C G Afo T(s) (x). 

Recalling (i) in Proposition 13. l\ if C G . (x) then h(x,C) — 0. Assume by contradiction that 

h(x,() < 0. Then, by using (ii) in Proposition 13. II there exists a > such that 

(a(,-l)EN epm (x,T{x)). 

Since A r epi ( T )(x, T(x)) is convex, we have, for any A G (0, 1), 

v x := A(C,0) + (1-A)«-1) 

= (AC + (1-A)<A-1) e N epi{T) (x,T(x)). 



This implies 



- I >-! ) £ N epi{T) {x,T{x)), 



1-A 



,A G 8T(x). 



1 - A 

By Theorem 5.1(b) in |26j . we have 

i.e., AC + (1 — X)aQ = A — 1, for all A G (0, 1). Letting A — > 1~ in the above equality, we obtain 
h(x, C) = 0, which is a contradiction. Thus h(x, () = 0. 

Conversely, let ( G N-j^.^Jx) be a nonzero vector such that h(x,() = 0. Applying Proposition 
13. H we see that (C,0) G N ep UT)(x,T(x)), which says exactly that ( G d°°T(x). The proof is 
concluded. 

□ 



4. Linear systems in 1 w 

This section is devoted to the study of non-Lipschitz points for the minimum time function to 
reach the origin for a linear autonomous system of the form: 

(4.1) x = Ax + Bu 

where A G M NxN and B G M NxM , 1 < M < N, and u G [-1, 1] M . 

We assume that (|4.ip is normal, i.e., for each column b{,i = 1, . . . ,M, of -B, the Kalman rank 
condition 

(4.2) rank[6i, Ab h A N - l hi\ = N 
holds. We set also 



(4.3) 

Of course, 1 < k < M. 



k = rankS. 
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We will first characterize the set S of non-Lipschitz points of T as 

S = jx G R N : there exist r > and £ G S^" 1 such that 

(4.4) M 

x = V / e A ^ r \ sign ((£, e At &i» dt and (C, &j) = Vi = 1, . . . , M \. 
i=l ^ J 

If = -AT, then 5 is empty. If < N, S is nonempty and we will prove also that S is (N — k)- 
rectifiable, with positive and locally finite % ^-measure. The positivity part will be contained in 
Section [JJ 

^From now on, we assume 

(4.5) k < N, 

Remark 4.1. It is easy to see that 

S = |x G R N : there exist r > and £ G S^" 1 suc/i that 
M „ r 

(4.6) X = Y, e ' M ^ si § n (K> e ~^» d *' 



i=l 



(G% r (x) and ((, e - Ar bi) = Vi = l,...,M}. 



We state first a technical lemma concerning an explicit computation of the minimized Hamiltonian. 
Before stating it, let us observe that condition f|4.2j) implies that the function 1 1— > (£, e~ At bi) is not 
identically for all i = 1, . . . , M. 

Lemma 4.2. Let r > 0, x G R N and ( G E> N ^ be such that 

M ,. r 

x = J2 e~ A % sign ((C, e- A %)) dt. 
i=i J o 

Then 

M 

(4.7) h(xX) = -^2\(C,e- Ar b l 



i=l 



Proof. We have 



h(x,Q = (C,Ax} + min (C,Bu) 

ue[-i,l] M 



M 

(C,M)+ min Y\C,Mi) 
Kl<i . r 



i=l,...,M 



1=1 



M „ r M 

C>E / A e - A *6 i sign((C,e-^))dt)-^|(C,& i )| 

1=1 ^° 8=1 

M , r 

J] / (C,Ae~ At b i siga{(C, e - At b i )))dt-\(C,b i )\ 
i=i ^° 

A/ 

^/tj(x,c). 



i=l 







U + t 
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Set gi(t) := (C,e~ At bi), t > 0. Then being not identically zero, gi vanishes at most finitely many 

times in [0, r], say at < t\ < ■ ■ ■ < t^ < r. We have, for i = 1, . . . , M, 
,. r 

hi(xX) = / -&(i)sigu(flj(i))dfc-|(C,&i>| 
Jo 

9i(t) sign(gi(t))dt 9i(t) sign 

3=1 Jt i 

- / ^(*)sign(ft(t))dt-|(C,6i)| 

= (ffi(0) - si g n {^9i yjj) + ~ 9i(tj+l)) sign 

+ (#(**) -0i(r)) sign f# (~^~~)) ~ KC.MI 
= 5i(0)sign - 0i(r)sign ^ f ~^y~~ J J ~ 5 r i(0)sign(^(0)). 

If 5i(0) / and ^(r) / 0, then sign(ft(0)) = sign [ gi (£■)) and signer)) = sign ( gi (^)). Thus 
hi(x,() = -\gi(r)\ = -\((,e~ Ar bi)\. Analogously, if ft(O) = and g^r) = 0, then hi(x,C) = = 
-|^(r)|. If gi(0) / and gi (r) = 0, then hi(x,Q = -\g%(r)\. Finally, if ^(0) = and gi (r) / 0, 
then h{(x,Q = — \gi(r)\. In all cases, we have 

h i (x,() = -\g i (r)\ = -\(C,e- Ar b t )\, 

and (jHZD follows. □ 

Remark 4.3. The characterization (|4.6p . thanks to Lemma \4.2\ and Theorem ] 2. 31 implies also that 

S = {x e R N : 3r > suc/i f/jai x G bdry7£ r and £ € S^ 1 n A^ r (x) /or w/iic/i /i(x,C) = 0}. 

The computation of the Hamiltonian contained in ()4,7p permits to prove the following charac- 
terization of non-Lispchitz points of T. We recall that, under the assumption (|4.5p . the set 5 is 
nonempty. 

Theorem 4.4. Let iel JV \{0}. Then T is non-Lipschitz at x if and only if x 6 5. Moreover, S 
is invariant for optimal trajectories of the reversed dynamics having vanishing Hamiltonian. 

Proof. Let i / be a non-Lipschitz point of T and set r = T(x) > 0. We recall that by Theorem 
3.7 in [13] epi(T) has positive reach. Therefore, by Propositions 13.31 and 13.41 there exists C £ 
§AT-l n jv^jc) suc h that C) = 0. 

Let u(-) be the optimal control steering x to the origin in the minimum time r. Then u(t) = 
u(r — t),0 < t < r, steers the origin to x in the optimal time r for the reversed dynamics x = 
—Ax — Bu,u G [—1, 1] M , namely 

M ,. r 

- A ^-%Ui(t)dt, 



I e- A ^Bu(t)dt = -f] f e 
Jo - =1 Jo 



where b{ are the columns of B and u = ({t 1} . . . , «m)- 
By the Maximum Principle, 



«i(t) = -sign ((C, e- A ^bi)) , i = 1, . . . , M. 
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Therefore, 



M rr 
M ,. r 

(4.8) = J2 e- A \s\gu{{le- A %))dt. 

„■ i Jo 



i=i 

Since h(x,() = 0, dUTJ) yields 

(C, 6-^=0, i = l,...,M, 

and the proof of x G S is concluded, recalling (|4.8p and (|4.6p . 

Conversely, let x G 5 and set r = T(x) > 0. Then by (|4.6p there exists £ G S^ -1 n Nfi r (x) such 
that (C, e" Ar 6i) = for all i = 1, . . . , M, and 

x = V / e-^&j sign ((C, e- Ai 6,)) dt. 
i=l ^ 

Recalling (|4.7p . h(x,Q) = 0. Therefore, by Propositions 13.31 and 13.41 T is non-Lipschitz at x. 

The last statement is an immediate consequence of Remark l4.1i In fact, the alternative espression 
of S given in (|4.6p . together with the Maximum Principle (see Theorem [23]) shows that every x € S 
is the endpoint of a time optimal trajectory for the reversed dynamics with vanishing Hamiltonian, 
starting from the origin. The proof is concluded. □ 

We prove now a rectifiability property for S, which is the main result of this section. 

Theorem 4.5. Let S be defined according to \4-4\ ) an d let k be given by \4-ty - Then S is closed and 
countably (N — k)-rectifiable. More precisely, for every r > there exist countably many Lipschitz 
functions fj : R N - k - 1 -> R N such that S n bdryft r C U j f j (R N - k - 1 ). 

Proof. The characterization of S contained in Remark 14.31 implies immediately its closedness. 
We now deal with the rectifiability of S. Observe that the set 

(4.9) Z := {C G S^ 1 : (C, h) = 0, Vi = 1, . . . , M} 
is a N — (k + 1) manifold. We define, for i = 1, . . . , M, 

gt(t,0 = ((,e A %), (ez,t>0, 

*<(r,C)= I" e A{t - r) b iS ign(g+(t,0)dt 
Jo 

and 

(4.10) Sj = {x G M N : there exist r > and ( £ Z such that x = $j(r, C)}- 
We claim now that 

(4.11) each S, is contained in a countable union of Lipschitz graphs of N — k variables. 

To this aim, we fix the index i and drop the corresponding subscript for the sake of simplicity. 
We set the following definitions. Fix r > 0. For every (N — l)-tuple of nonnegative integers, 
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we define 

Zj = G Z : g + (t,() has in the interval [0, r] exactly 
ji zeros of multiplicity 1, 



JN-1 



zeros of multiplicity N — 1 1 . 



We set also |j| = ji + • • • + 3n-i and observe that, thanks to (|4.9|) . we can consider only j's such 
that |j| > 1. Moreover, for any positive integer d and j G N^ -1 with |j| > 1 we define 

Zf = jce Zj : min{|ri - r 2 | : <? + (ti,C) = 5 + (r 2 ,C) = O.n / r 2 } > 1 



Invoking Lemma 3.2 in [To], we obtain that 



IJI=1 




We define finally the map 
(4.12) 

and, for all j 6 N^ -1 , the sets 



Y : Z ->• L^O.t) 

C H> sign(5+(-,C)), 



2f ± = {<eZf: ( M sign( 5 + (t,C)) = ±l}. 



We fix now j € N w . If |j| = 1, then K(C)(t) = ±1 for all ( € Z h t G (0,r], and so K is locally 
Lipschitz in Zj. We claim that Y is locally Lipschitz also in Z. ' + and in Z. '~ for each j G N^ -1 . 
The argument for Z- ,+ and Z. ' is the same, so we perform it only for Z- ,+ , |j| > 1, d > 1. 

So, fix |j| > 1, d > 1 and Co £ ^/' + - Let • • • , t\j\ be the zeros of g + {-,(o) m [0, r], each one 
with multiplicity rrih, h= 1, . . . , |j|. 

By continuity and the implicit function theorem, for each h = l,...,|j| there exist a compact 
neighborhood Vh of Co, a neighborhood Ih of th and a C 1 -function ip^ : Vh — > Ih such that 

(4.13) __ 5 +(t, c) ^o V(t,()e4xi4 

and 

(4.14) | (C, t) &Vh x Ih '■ ^ ZT g + (t, C) = j = graph(^). 

The neighborhoods can be taken disjoint and satisfying |J^| < We choose now 1/ = V(Co) Q 
njj^jVh with the further requirement that for all ( <E V, the set {t G [0,r] : g + {t,() = 0} is 
contained in Uj|Li ^h- Since |J^| < the function g + (i, C) has at most one zero in each Ih- 

Set Vj(Co) = V n Without loss of generality, we may assume that Z?' + is contained in a 

finite union of such Vj(-)> say = U^^j(O)- We write the functions corresponding to Vj((^) as 
</>h(C)> h = 1, . . . , |j|, and observe that each ^ is Lipschitz continuous on Vj(Q), say with Lipschitz 
constant L h . We denote also the intervals corresponding to V$(Q) as ifa, h = 1, . . . , |j|. Of course, 
some of the Vj((^)'s may be the singleton {Q}, and in this case everything trivializes. Fix now an 
index I. 
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We claim that, for each £ £ Vj(Q), g + (-,() has a zero of multiplicity rrih exactly at <p e h ((), 
h = 1, . . . , |j|, and does not have other zeros in [0, r]. Indeed, by construction for each £ S Vj(Q) 
all zeros of C) are contained in (J/!Li ^I- Let k be the largest index k such that 7^ 0. Again 
by construction, for each Q G Vj(0) the map t h-> <7 + (t,C) has exactly j K zeros of multiplicity re. 
Moreover, such j K zeros must belong to the same intervals ii to which the j K zeros of multiplicity n 
°f 9 + ('iCi) belong, since in all other intervals we have at least one nonvanishing derivative of order 
< k — 1. Owing to (|4.14p with = k, such zeros must occur at <^(C)> f° r the corresponding 
index /i. Let now k\ be the largest positive integer < k such that j Kl > 0. By definition of Vj(Q), 
for each £ £ ^j'(C^) the map t \-t g + (-, Q does not have zeros of order k, with k\ < k < n and must 
have exactly j K1 > zeros of multiplicity K\. Such zeros cannot belong to the intervals to which the 
re-zeros of <? + (-,C) belong, since such intervals already contain a zero; on the other hand, by (|4.13|) 
they must belong to the same intervals il to which the zeros of multiplicity k\ of g + (-,Ce) belong, 
and therefore they must occur at 9?^(C); f° r the corresponding index h. An analogous argument 
can be performed for all further indexes k < k\ such that jk 7^ 0. Therefore the claim is proved. 

We are now ready to show that Y is Lipschitz on Vj(C^)- Indeed, fix the index t and let Ci> C2 £ 
V^Q). Then 

Li I 

unc2)-nci)LW) < 2 kto) - <4(coi 

h = 1 

rrih is odd 

<2 y, 41K2-C1II, 

h = 1 
rrih is odd 

which proves the claim. 

The Lipschitz continuity of Y on each V;(Q) implies immediately that, for all fixed r G [0, r], the 
function £ i-> $(r, C) is Lipschitz in the same set. On the other hand, the function r 1— > <&(r, Q is 
immediately seen to be Lipschitz on [0, t}. Consequently, the set £ defined in (|4.10p is contained 
in a countable union of Lipschitz graphs of iV — k variables. The (N — A;)-rectifiability of S now 
follows easily and the proof is concluded. □ 



5. Nonlinear systems in R 2 
This section is devoted to the study of non-Lipschitz points of T for the nonlinear system 



(5.1) 



x(t) = F(x(t)) + G(x(t))u(t), 

u(t)e[-l,l] M , 

s x(0) = x, 

where the state x is in M 2 and M is either 1 or 2. 
The assumptions are the following: 

1) F : M. 2 — > M 2 and G : M 2 — ?• M.2xM are of class C 1,1 and all partial derivatives are Lipschitz 
with constant L; 

2) F(0) = 0; 
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3) Tank[Gi(0),DF(0)Gi(p)] = 2 for i = 1, . . . , M, where we mean G = Gi if M = 1 and 
G= (Gi,G 2 ) if M = 2; 

4) ,00(0) = 0. 

Theorems 5.1, 6.2 and 6.5 in [15] yield that there exists T > 0, depending only on L, DF(0), and 
G(0), such that for all < r < T, 

a) 7?. r is strictly convex and for all x G bdry7£ r there exists a unique optimal control u(-) 
steering x to the origin in the minimum time r, and u(-) is bang-bang with finitely many 
switchings, 

b) every x G 7£ T is optimal (the definition of optimal point was recalled in Section ETTj) , 

c) epi(T) has locally positive reach. 

We prove here a result which is the nonlinear two dimensional analogue of Theorem 14.51 Fix 
< r < T and define 

(5.2) S = jx G 1Z T : 3( G S 1 n Afa T(a!) (a) such that h(x, C) = o} . 

Recalling Propositions 13.31 and [3T4"1 S is exactly the set of non-Lipschitz points of T within 1Z T . 
We show first that S is invariant for a class of optimal trajectories and then that it is countably 
'W 1 -rectifiable. 

Proposition 5.1. Let S be defined according to (|5.2p and let F, G satisfy the assumptions 1) - 
4)- Then S is invariant for optimal trajectories. 

Proof. We wish to prove that if x G S and x(-) is the optimal trajectory steering x to the origin in 
the minimum time T{x) then x(t) G S for all < t < T(x). In fact, let u(-) be the corresponding 
optimal control and set u(t) = u{T(x) — t),0 < t < T(x). Let x(-) be the solution of the system 



(5.3) 



±{t) = -F(x(t)) - G(x(t))u(t), 
x(0) = 



and let (sS'n N nT (x) be such that h(x, C) = 0. 

Claim. The solution \(t), t G [0, T{x)\ of the adjoint system 



(5.4) 



\{t) = \(t)(DF(x{t)) + DG(5:{t))u{t)) 
\(T{x)) = C, 



satisfies the following properties: 

(i) Ui(t) = sign((A(i),-G i (x(t)))) fora.e. t G [0, T(x)], i = 1, . . . , M, 

(ii) = h(x(t),\(t)) = (F(x(t)),\(t)) ~ E£i \{Gi(x(t)),\(t))\ for all t G [0,T(x)], 

(iii) ^ A(t) G ^ T( , (t)) (x(t)), for all i G [0,T(s)]. 

Proof of the Claim. We recall that under our assumptions TLt{x) is strictly convex. In particular, 
Nji t ^Jx) is the convex hull of its exposed rays (see [231 P-163] and [231 Corollary 18.7.1, p. 169]). 
Therefore let ( ^ belonging to an exposed ray of Nn T , s Jx). Recalling (b) in Proposition 13.11 

there exists a < such that , ££jS2 belongs to an exposed ray of AL Di m(x, T(x)). 

VllCII +o- v ' 

By Theorem 4.9 in [12], there exists a sequence {x n } C dom(LT) such that and 

lim ( jDT ( x «)'- 1 ) _ (M 



y/\\DT(x n )\\* + l vWT 



0~ 



NON-LIPSCHITZ POINTS AND THE SBV REGULARITY OF THE MINIMUM TIME FUNCTION 



17 



Let u n = (ux >n , . . . ,UM,n) be the optimal control steering the origin to x n . Since N-r, t ,. (x n ) is the 
half ray M + DT(x n ), for n large enough, then Pontryagin's Maximum Principle yields that 

(5.5) Ui,n{t) = sign ((X n (t), —Gi(x n (t))}) , a.e. t G [0,T(z)], i = l,...,M, 
where x n (-) is the solution of 

| y = - G{y)u n 

{ y(o) = o, 

and A n is the solution of 

f A(t) = A(t)(L»F(x n (t)) + J DG(x n (t))u n (t)), a.e. 
\ X(T(x n )) = Cn G R+£>T(x n ). 

Since all controls u n are bang-bang with a finite number of switchings independent of n, up to a 
subsequence we can assume that u n (-) (where we have put u n (t) = for t G (0, T(x) — T(x n )) if 
T(x n ) < T(x)) converges pointwise a.e. to some admissible u : [0,T(x)] -)■ [-1,1] M . Let x (-) 
be the solution of (|5.3p with no in place of u. Since obviously xo(T(x)) = x, by the uniqueness of 
the optimal control we have that uo(t) = u(t) a.e. on [0, T(x)]. Up to another subsequence, we 
can assume that x n (-) converges uniformly to x(-) on [0, T(x)], and A n (-) converges uniformly to 
A(-) on [0, T(x)]. Then A(-) is the solution of (|5.4p with £ in place of Recalling (|5.5p . the above 
convergence properties imply that 

(5.6) Ui(t) = sign ({\(t),-Gi(x(t)))), a.e. t€ [0,T(x)\, i = l,...,M. 

Let now Ci) C2 G S 1 belong to exposed rays of N-ji T , a , (x) and let a, ft > be such that C = a^i +(3(2- 
Let Ai(-) (resp., A2(-)) be the solutions of (j5.4[) with (1 (resp., (2) in place of (. By (|5.6p . we have, 
for a.e. t G [0,T(x)], that 

Ui(t) = sign ((Ai(t), -G<(£(t)))) = sign ((A 2 (t), -G?i(5(t)))) , i = 1, . . . , M. 

Therefore, for a.e t G [Q,T(x)], 

Ui(t) = sign((a~Ai(t) + p\ 2 (t),-Gi(x(t)))^ 

= sign((A(t), -Gi(x(t)))), i = l,...,M, 

which proves (i). 

To prove (ii), observe that the fact that h(x(t), A(f)) is constant follows in a standard way from the 
maximization property (i) (see, e.g., Corollary 6.4 in [IS]). Since h(x(T(x)), X(T(x))) = h(x, Q = 0, 
(ii) is proved. 

Statement (hi) again follows from the maximization property (i) (see, e.g., Remark 5.2 in |15|). and 
the proof of the Claim is concluded. 

We now complete the proof that S is invariant for optimal trajectories. To this aim, fix x G S, 
together with £ G S 1 D N-ji t ^ such that h(x,() = 0. By the above claim, the never vanishing 
adjoint vector A(-) which is the solution of (|5.4p is such that h(x(t), A(t)) = and X(t) G Nn t (x(t)) 
for all t G [0, T(x)], which shows that each point x(t) of the optimal trajectory x(-) steering the 
origin to x belongs to S. The prove of the invariance of S is complete. □ 

Theorem 5.2. Under the same assumptions of Proposition HOI the set S is countably Ti 1 -rectifiable. 
Moreover, for all x G S there exists 5 > such that 

(5.7) tf(SnB(x,(5)) > 0. 
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(5.8) 



Proof. In order to prove the rectifiability property of S, it is enough to show that, if S is nonempty, 
then it consists exactly of two optimal trajectories of the reversed dynamics 

f x(t) = -F(x(t))-G(x(t))u(t),ue[-l,l} M ,te[0,r], 
\ s(0) = 0. 

Let i£S together with C G S 1 n N-j^.^ be such that h(x, () = 0. Let u(-) be the optimal control 
steering the origin to x and let x(-) (resp., A(-)) be the corresponding optimal trajectory (resp., 
adjoint vector, the solution of (|5.4|) ). Set Co = A(0) / 0. 

We assume now that M = 1, i.e., the control is scalar. Since the Hamiltonian is constant along 
the optimal trajectory x, we have that 

|(G(0),Co>| = fc(0,Co) = 0(= MO, -Co))- 

We now prove that each one of the vectors Co an d —Co determines uniquely an optimal trajectory 
of (|5.8p contained in S. In fact, for every optimal trajectory x(-) of (|5.8p . with a corresponding 
adjoint vector A(-), we can define the switching function 



gtx(t) = (-G(x(t)),X(t)). 
/,. a \vm • ' • • su/ = and g+ x ( 



Of course, 3+ A (0) = (-G(0),±Co> = and 5+ A (0) = ^(DF(O)G(O), Co)- The last expression is 



nonzero, due to the assumption 3), so that in a neighborhood of t = 0, the sign of <?Z\(') is 
uniquely determined by ±Co- Therefore, in a neighborhood of t = the optimal control is uniquely 
determined by sign(<jl\ (•)), by the Maximum Principle, and so there are exactly two optimal 
trajectories of (|5.8p which belong to 5 in a neighborhood of t = 0. Since at every zero of <?Z\ (•) the 
derivative g~\(') is nonvanishing (see, [151 Sections 3.2 and 5]), the optimal control can be uniquely 
extended up to the time t = r. The proof is now complete for the case of a single input. 

To conclude the proof of the rectifiability, let M = 2. The condition h(0, Co) = means that the 
system of equations 

f <Gi(0),Co> =0, 
\ (G 2 (0),Co> =0. 

has nontrivial solutions. So, if Gi(0) and Cr2(0) are linearly independent, then S is empty. Oth- 
erwise, both components of the optimal controls are uniquely determined by the sign of the corre- 
sponding switching functions, exactly as for the single input case. 

The propagation property ()5.7[) is an immediate consequence of Proposition 15.11 The proof is 
concluded. □ 

6. The SBV regularity of T 

As a consequence of the results contained in Sections 4 and 5 we prove the SBV regularity of 
the minimum time function T. We recall first some properties of functions with bounded variation, 
and next we collect some known results on functions having epigraph with positive reach. As it 
was proved in [13] and in [15] . the minimum time function has this property under the assumptions 
taken in Section 4 or in Section 5. 

Let f2 C be open. We say that a function / G L^ oc (Q) has locally bounded variation, and 
we denote this fact by / G BV\ oc (£l), if for every ball A C ft the distributional derivative of / 
in A is a finite Radon measure (see, e.g., [U Definition 3.1]), which we denote by Df. We write 
Df = D a f + D s f, where D a f is absolutely continuous with respect to Lebesgue measure, and 
D s f is singular. The singular part D s f can also be decomposed into the jump part, D 3 f, and the 
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Cantor part, D°f (see, [H Section 3.9]). In the case where / is continuous, like in the case / = T 
under our assumptions, the jump part obviously vanishes. 

Definition 6.1. (see, e.g., [1, Section 4.1],) We say that f G BV\ oc (Sl) is a special function of 
locally bounded variation, f G SBV\ oc (Sl), if the Cantor part of its derivative D c f vanishes. 

It is our aim, in this section, to prove that under the assumptions of Section 4 and 5, the Cantor 
part D C T vanishes, and so T is a special function of locally bounded variation. 
We state some further results. 

Proposition 6.2. (see [Tj Proposition 4.2]) Let f G BV(Sl). Then f G SBV (SI) if and only if D s f 
is concentrated on a Borel set a -finite with respect to % N ~ l , in particular, if it vanishes outside a 
countably % N ~ l -rectifiable set. 

Recalling the definition of non-Lipschitz points given in Section 3 (see Definition I3.2|) . we obtain 
the following result. The notation /ing means the restriction of the measure [i to the set E. 

Proposition 6.3. Let SI C R N be open and let f G SBV\ oc (Sl). Let 

K = {x G SI : f is non-Lipschitz at x}. 

Then D s f mK = 0. 

Proof. By definition, / is locally Lipschitz in the open set SI \ K. Therefore Df is absolutely 
continuous with respect to Lebesgue measure C N in SI \ K (see, e.g., [TJ Proposition 2.13]), i.e., 
D s f L n\K = 0. □ 

Consequently, by putting together the two previous Propositions, we obtain 

Corollary 6.4. Let SI G 1^ be open and let f G BV\ oc (Sl). Assume that the set of non-Lipschitz 
points of f be countably Ti N ~ l -rectifiable. Then f G SBV\ oc (Sl). 

We are now ready for the main results of this section. 



Theorem 6.5. Consider the linear control system fl^.i[ ) under the assumption M-ty . Then the 
minimum time function T to reach the origin satisfies T G S-BVioc^^)- 

Corollary 6.6. Let 1 < J < N be the smallest integer such that 



B,AB,...,A 



J-i 



N. 



rk 

Then for every 1 < p < N(N — J + l)/(A r — J) we have 

TGW^(R N ). 

Theorem 6.7. Consider the nonlinear system \5. 1\) under the assumptions 1) - 4) stated in Section 
5. Then there exists T > depending only on G(0), DF(0), and on the Lipschitz constant L of 
DF and DG, such that T G S BVi oc {mt(R r )) ■ 

Proof of Theorem \6.5\ and \6. 7[ The statements follow immediately by putting together Theorem 
2.21 and Theorem 14.51 (resp.. Theorem I5.2p and Corollary 16.41 □ 

Proof of Corollary 1 6. (A Since T is continuous and belongs to S^Vio^l^) then its distributional 
derivative DT is a locally summable function. Moreover, it is well known that T is Holder continuous 
with exponent 1/J (see, e.g., [2J Theorem IV. 1.9] and references therein). The statement then 
follows by applying standard results on Sobolev spaces (see, e.g., Theorem 3, p. 277, in [I?]). □ 



20 



GIOVANNI COLOMBO, KHAI T. NGUYEN, AND LUONG V. NGUYEN 



7. Propagation of non-Lipschitz singularities 

This section deals with a lower estimate of the dimension of S for the linear case. We show 
that the 'H Ar_fc -rectifiability of S is indeed optimal, in the sense of Theorems 17,11 and 17.21 below, 
at least for a small time. Those statements can be seen as propagation results for singularities of 
non-Lipschitz type for the minimum time function. 

We consider the linear system (|4.ip under the assumptions (|4.2|) and (|4.3p . Let N > 2 and define, 
for r > 0, <S(t) = S n bdry7£ T . We assume that k < N — 1, otherwise 5 = 0. 

Theorem 7.1. There exists f > 0, depending only on A, B, N, satisfying the following the property: 
for all t < f and all Tl N ~ 2 -a.e. x S 5(r) there exists a neighborhood V of x such that 

(7.1) ^"(Fn5(r))>0. 

Proof. We divide the proof into some steps. 

We will use a result which was proved, e.g., in [2] (see the proof of Lemma 8). The statement is as 
follows. 

There exists r > 0, depending only on A, B, N, such that for every C G 8 N ~ 1 the switching 

(7-2) 

function g(-, () = (£, e b) has at most N — 1 zeros in [s, s + r] for every s > 0. 

Claim 1. The statement of the Theorem holds true in the case B = b, a vector. 

Proof of Claim 1. Let f be given by (|7.2p . Fix < r < f and let < j < N - 2. We say that 

x £ S(t) belongs to 5? (t) if there exist times < si < S2 < ■ • • < sj < r such that 

x = ± f e~ As b-/(s)ds, 
Jo 

where 



7(s) 



1 if < s < si, 

— 1 if Si < S < S2, 



(-IP if Sj <S<T. 

In other words, the optimal control steering the origin to x for the reversed dynamics has exactly 
j switchings in the interval (0,r). 

Step 1. There exists f > 0, depending only on A, b and N, such that if < si < S2 < • • • < Sj < f 
then 

(7.3) rank [e~ Asi b, e~ As2 b, ■ ■ ■ , e~ As i b] = j. 

In order to prove ([731), for s > and set g(s) = ((, e~ As b) and H = {( G R N : g( Si , () = 

0, i = 1, . . . ,j}, < si < S2 < • • • < Sj < t. We claim that dimH = N — j. Indeed, if g(si, () = 0, 
then 

= (C,b-A Sl b)+o(r), 

so that (£, b) = since r can be chosen small enough. Furthermore, if j > 1, there exists s\ S (s\, S2) 
such that J^g(s~i, £) = 0, which in turn implies 

= -{(,Ab-A 2 s 1 b) + o(r), 
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so that (C,Ab) = since r can be chosen small enough. The same argument provides times Sj, 
i = 2, . . . j — 1 such that 

= ^(s i ,0 = (C,A i b) + O(r), 
i.e., {(,,A l b) = 0. The proof is completed by invoking the rank condition (|4.2p . 

Step 2. If < j < iV — 2, then for all < r < f the set 5j(r) is the union of two smooth 
parametrized j-surfaces. Actually we are going to prove that {x G Sj(t) : x = C e~ As bj(s)ds} is 
a smooth parametrized j-surface, the other case being entirely analogous. 
Indeed, we have 

x = r e~ As bds + ^(-1)* e~ As bds + (-iy f T e~ As bds, 

JO i=1 J Si Jsj 

where < si < . . . < sj < r. Observe that §§- = 2(-l) i+1 e - As 'b, 1 < i < j, and by Q the 

matrix ( Mf- ) has rank j in the open set {(si, . . . , Sj) E (0, r)- 7 : < . . . < Sj}. The proof of 

V 1 / i=i,— ,j 

Step 2 is concluded. 

Set now 5o(t) = {± f Q T e~ As bds}. By the Maximum Principle, owing to (j7.2|) we have that 

AT-2 

5(r) = |J S 3 (r) 

for all < r < f and the union is disjoint. In particular, Step 2 implies that for all such r 
(7.4) n N - 2 {S(r)\S N ^(r))=0 

and that (|7.ip holds at every point x G 5at_2(t). The proof of Claim 1 is concluded. 
Claim 2. The statement of Theorem 17.11 holds in the general case. 

Proof of Claim 2. Let < r < f be given and fix x G S n bdry7?. r , together with the optimal 
control u = (ui, . . . , «m) steering the origin to x in time r by the reversed dynamics. Assume that 
Ui has exactly Ki + 1 zeros, < Ki < N — 2, at times 

= 4 < s\ < . . . < 4. < r. 

Then, recalling (|4.3|) . we have 

k < rank{6i, e As h h . . . , e As ^6i : 1 < i < M) < N - 1. 
For j = 0, . . . , N — (1 + k), let Sj + k(r) be the set of all x G <S(r) such that 

rank{e As i&j, . . . , e As ^ 6j : 1 < i < M} = j + 1. 

Observe that 

N-(i+k) 
S(t)= [J S J+k (r) 

3=0 

and the union is disjoint. Moreover, by arguing exactly as in Steps 1 and 2 in the proof of Claim 
1 above, we can see that each Sj + k( T ) 1S a union of finitely many disjoint smooth parametrized 
^'-surfaces. Thus, 

n^-^Sj+kir)) =0 Vj = Q,...,N-(2 + k) 

and 
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Therefore, for H N ( fc+1 )-a.e. x G <S(r) there exists a neighborhood V = V(x) such that 

n N -( k+1 \vns(r)) >o. 

The proof is now complete. □ 

By combining the above result with the invariance statement contained in Theorem l4.4l we obtain 
immediately the following 

Theorem 7.2. Let f be given as in Theorem \7.1\ Then for 7i k -a.e. x G S such that T(x) < f 
there exists a neighborhood V of x such that 

H N - k (vnS) > o. 

Proof. Fix < r < f and let E be a subset of S n bdry7£ T with full "H^^'+^-measure with the 
property f)7. 1 1) . Then the optimal trajectories for the reversed dynamics through each point of E 
from a subset of S with full "H^^-measure. The proof is complete. □ 

Remark. The statement of theorems 17.11 and 17.21 are somewhat unnatural for linear systems, as 
they are valid only for small times. The proof for arbitrarily large times requires an analysis of 
higher order and of linearly dependent zeros of the switching function, which we are not yet able 
to conclude. 

Observe finally that in the nonlinear two dimensional case the propagation result is contained in 

(S3J). 
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